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We show that the usefulness of the thermal state of a specific spin- lattice model for measurement- 
based quantum computing exhibits a transition between two distinct "phases" - one in which every 
state is a universal resource for quantum computation, and another in which any local measurement 
sequence can be simulated efflciently on a classical computer. Remarkably, this transition in compu- 
tational power does not coincide with any phase transition, classical or quantum, in the underlying 
spin-lattice model. 
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I. INTRODUCTION 

Determining whether properties of the ground or ther- 
mal state of a coupled quantum many-body system can 
be efficiently simulated on a computer is a central ques- 
tion in many-body physics. Consider a spin lattice model 
that is unitarily equivalent to a lattice of uncoupled spins; 
in such a model, the spectrum and the behaviour of 
simple correlation functions can be solved analytically. 
However, such simplicity can belie an underlying com- 
plexity, as the calculation of more general observables 
on the thermal state of this system can be computation- 
ally intractable. For example, one might want to sim- 
ulate the outcomes of a sequence of local (single-spin) 
measurements on a subset of spins; such measurements 
are used to reveal "hidden" order such as long-range en- 
tanglement [l|]. Although numerical techniques such as 
Monte Carlo methods [l| may be efficient in some in- 
stances, one can also devise models for which the out- 
comes of such measurements performs a quantum com- 
putation [3, |3| for which no classical simulation algorithm 
is expected to be efficient. Little is known about how to 
characterize the classical simulation complexity of such 
general observables, even in the simplest models. 

Here, we provide a remarkable example of a solvable 
model wherein the efficiency of simulating general observ- 
ables of the thermal state undergoes a transition, separat- 
ing the model's parameter space into two distinct regions. 
In one region, at high temperatures, we present an ex- 
plicit algorithm that can efficiently simulate the outcome 
of any local adaptive measurement sequence. In the other 
region, we prove (based on standard complexity assump- 
tions^) that no such efficient simulation algorithm exists. 



Remarkably, there is no phase transition in this model, 
classical or quantum, that could serve to demarcate these 
two regions. 

Specifically, we make use of the cluster state model of 
measurement-based quantum computation (MBQC) [3, 
y], wherein a fixed resource state is subjected to local 
measurements [1]. We consider a spin-lattice model for 
which the unperturbed ground state is a cluster state 
but which is subject to both thermal noise and a uni- 
form external field. We introduce a method for efficient 
classical simulation of local measurements on sufficiently 
high-temperature equilibrium states of this model. Our 
method significantly generalizes the separability result 
for bipartite partitions of thermal cluster states (5[ by 
combining ideas from mixed-state entanglement and per- 
colation theory. For suSiciently low temperatures and 
small external fields, we describe a method for obtain- 
ing cluster states on which MBQC can be performed by 
making use of the error-correcting thresholds of p, |7| . 

Our paper is structured as follows. In Sec. |IT1 we 
present our completely-solvable model. We provide a 
lower bound on the temperature for which MBQC on 
the thermal state can be efficiently classically simulated 
in Sec. Illli and an upper bound on the temperature for 
which the thermal state is a universal resource for MBQC 
in Sec. IIVI Conclusions are presented in Sec. |Vl 



^ It is believed that local measurements on such states are impos- 
sible to efficiently simulate on a classical computer, since, if this 
were possible, any polynomial time quantum algorithm could be 
performed efficiently on a classical device. 



II. CLUSTER-STATE MODEL 

The model we consider is the cluster state Hamilto- 
nian [5| with the addition of a local magnetic field in the 
z-direction. A cluster state is a highly entangled state of 
two- level systems on a lattice £ @, 3 • It can be charac- 
terized as the unique +1 eigenstate of a set of commuting 
operators Ki — XiY\-^Zj, where Xi (Zi) is the Pauli X 
(Z) operator at site i and where j~z denotes that j is 
connected to i by a bond in the lattice £. The cluster 
state Hamiltonian [5|, |8| is a model of interacting spins 



with H = —^ SiGz; -^« ^'-'^ which the cluster state is the 
unique ground state, and which possesses an energy gap 
of A between the ground and first excited states. We will 
consider a modified Hamiltonian with a local Z field at 
each site: 



He = -#E.e£(cos K, + sin 9 Zi^ 



(1) 



for < < 7r/2. With this parameterization, the ground 
state energy gap of the Hamiltonian is A for all values 
of 6. The parameter 9 quantifies the relative strength of 
the local magnetic field term. For = we recover the 
cluster state Hamiltonian, while for 9 = tt/2 the spins 
are uncoupled and the ground state is a product state. 

The thermal state of this model can be found straight- 
forwardly for any strength of external field and temper- 
ature T because the system is unitarily equivalent to a 
lattice of uncoupled spins Q, as follows. The controUed- 
phase gate CZ = diag(l, 1, 1, —1) can be applied to any 
bond in the lattice; because of its symmetry, it does not 
depend on which qubit is the control and which is the 
target. Also, as the gate is diagonal in the Z-basis, its 
action on different bonds all commute with each other. 
Therefore, we can define the unitary CZ^ which is given 
by the product of controUed-phase gates on all bonds of 
the lattice. This unitary operator decouples the system, 



t\- 



iCZc)Z,{CZl) ^ Z, , 
{CZc)K,{CZ[)^X,, 



(2) 
(3) 



for all i G C. We can therefore express the transformed 
Hamiltonian as 



{CZc)HgiCZl 



fE,;e£(coseX, 



sin 9 Zi 



(4) 



i.e., as a sum over single spin observablcs. Note that 
CZ£ is a non-local unitary that creates entanglement, so 
that the ground state of He can be highly entangled even 
though Hq is unitarily equivalent to a model of uncou- 
pled spins. Note that this solution does not imply that 
MBQC can be efficiently simulated for any value of 9 or 
T, because the statistics of measurement outcomes may 
still be hard to calculate. 

This solution yields an explicit expression for the ther- 
mal states of ([1]). Define the single spin states 



= cos(6l/2)|+)+sin(6'/2)|-) 



(5) 



with |±) the eigenstates of X. The ground state of ([T|) 
is obtained by applying CZ£ to the product state \9)®'~ 
where each spin on L is prepared in the state \9). We 
recognize the ground state of our model as one of the 
states whose usefulness for MBQC was considered in jlO| . 
The thermal state p(/3, 9) is obtained by applying CZ£ 
to the thermal states 

\{I + tanh(/3A/2)e'''^'/2x,e-*''^-/2) ^ ^g^ 

of each independent spin, where /3 — (fc^T)"^. 



The model is unitarily equivalent to uncoupled spins, 
and so there are no phase transitions, quantum or clas- 
sical. Specifically, in the thermodynamic limit, the free 
energy per site is analytic for the full range of T and 9. 
Despite this, we now demonstrate that this system ex- 
hibits a transition in its usefulness as a resource for uni- 
versal MBQC. First, we demonstrate that above a certain 
finite temperature it is possible to efhcicntly classically 
simulate any attempt at MBQC. We label this region of 
the parameter space of the model 'Region C". Second, we 
demonstrate that for 9 and T below some specific thresh- 
olds, the thermal state is a universal resource for MBQC; 
this is denoted 'Region Q\ 



III. REGION C: EFFICIENTLY SIMULATABLE 

We define classical simulation of an attempted MBQC 
as in Ref . [ll| : Let N be the number of two- level systems 
in a particular finite sized lattice. An MQBC on these 
qubits can be efficiently classically simulated if there ex- 
ists a randomized algorithm that can sample results of 
arbitrary single-qubit measurements (together with feed- 
forward) from the correct (quantum-mechanical) proba- 
bility distribution using resources that scale polynomially 
in N on a classical computer. 



A. Simulation of thermal states 

Certain Monte Carlo algorithms simulate properties 
of a thermal density matrix p(/?, 9) by sampling energy 
eigenstates with probabilities given by the Gibbs distri- 
bution. Here we will develop a stochastic simulation that 
efficiently samples from a different set of states |Cq)(Cq,| 
with probabilities Pq, such that 



p(/3,0) = ^P„|C„)(C„| 



(7) 



(We note that the states |Cq,) need not be orthogonal, 
and do not in general form a basis.) The states |Cq) will 
be be chosen such that, for sufficiently high temperatures, 
they have such small amounts of entanglement that they 
can be stored and manipulated efficiently on a classical 
computer. We conclude that it is not possible to perform 
MBQC on the states p(/?, S) in the regime of temperature 
in which this simulation succeeds. 

The first step in our simulation algorithm is to repre- 
sent the state p(/3, 9) as a 'projected entangled pair state 
(PEPS) [141; a class of states that generalizes the well- 
known valence bond solids [Tj|. In a PEPS representa- 
tion, as indicated in Figurema) for a square lattice, a set 
of d{i) virtual quantum systems (we restrict ourselves to 
considering virtual qubits) is associated with each site i 
of the lattice, where d{i) is the number of bonds emanat- 
ing from site i. For every pair of sites i and j connected 
by a bond, we can identify two virtual qubits (one at site 



i and one at site j) which we associate to this bond. A 
PEPS on the original lattice is obtained by placing every 
such bond in some state pbond, and at every site i some 
operator Ai maps the d{i) virtual systems onto the phys- 
ical system. (The operators Ai are isometrics, but are 
commonly referred to as projectors.) 

For the cluster state (and generally for PEPS rep- 
resentations), the pair of virtual qubits associated to 
each bond are assigned the maximally- entangled state 
Pbond = 1^2) (C2 1, with 



(a) 



\C2) - 4. 



V2 



|l>|-», 



(8) 



where |±) 



1 

x/2 



(|0) ± |1)). With this convention, the 



cluster state has a simple PEPS representation [l^l cor- 
responding to the choice of isomctry 



A=|0)(00...0| + |1)(11...1|, 



(9) 



at each site, with d{i) zeros (ones) in (00 . . . 0| ((11 . . . 1|), 
and the states |0) and |1) forming a basis for the resulting 
qubit at each site. 

In our simulation algorithm, we follow |5[ and choose 
more general states for the pairs of virtual qubits on a 
bond, allowing us to describe thermal states of our cluster 
Hamiltonian, both with and without the local field term. 
The zero-field thermal state p{(3, 0) is obtained as a PEPS 
if we choose 

/Obond(/3, = 0) = i(/ + UJeX ® Z){I + UJeZ ® X) , (10) 

on every bond, and 



yl==|0)(00...0| + |l)(ll...l|. 



(11) 



at every site [5|. For uje ~ I, the bond states are pure 
and maximally entangled, and the resulting PEPS is the 
cluster state. For uje = 0, the bond states are maximally 
mixed, corresponding to infinite temperature. To obtain 
p{P, 0), the parameters uje are chosen such that 



J|we = tanh(/3A/2) 



(12) 



where the product is taken over the bonds emanating 
from a particular site; for simplicity, we choose all Ue to 
be equal. The case of non-zero 9 is handled by a slightly 
more general ansatz for the bond states, as follows 

Pbond(/3,6i) = \{I+aX®Z+"/Z®I){I+aZ®X+-fI®Z) . 

(13) 
The details of how to choose a and 7 to reproduce p(/3, 6) 
is described in Appendix IB] 

The second step of the algorithm is to decompose the 
states pbond as an ensemble of states, such that the ele- 
ments of the ensemble are product states with high prob- 
ability. Every two-qubit state has a decomposition of the 
form [11 



Pbond =Pe|V')('0| +EuPs,mI'^1,A')(<^1,mI '^ I'^2,a<) ('/': 
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FIG. 1: (a) A PEPS representation of the thermal state on 
a square lattice. Bonds between virtual qubits, denoted by 
dashed lines, are in a mixed state, (b) An instance of this 
thermal state. Solid lines denote a maximally-entangled state 
on the bond, whereas no line denotes a separable state. 



where \^) is a pure entangled state, |0i,^) ® \4'2.,^J) is a 
pure product state, p^, Ps.fj. are probabilities, and the 
probability of getting a product state X^/^Ps.^ = 1 ~Pe is 
maximal over all possible decompositions. Following [l^ , 
we find that for pbond(/9, 0) the state \tp) is simply the 
two-qubit cluster state, and that 



{luI + 2lo,-1)/2. 



(15) 



(14) 



As expected, high temperatures yield small values for 
Pe and an ensemble that is largely made up of product 
states. For the case of general 9 this decomposition is 
readily performed numerically. 

Our key insight for simulating MBQC on a thermal 
PEPS state is that, armed with a decomposition such as 
Eq. (|14p. one can efficiently generate instances of an en- 
semble that gives the state p{f3, 9). If the temperature is 
sufficiently high, the members of this ensemble can be ef- 
ficiently stored and manipulated in a classical computer. 
We outline the simulation algorithm here; full details are 
given in Appendix \X[ 

On a given run of the simulation, the state for each 
bond is chosen to be either the pure two-qubit entan- 
gled state IV') or one of the product states \4>i,i_i) ® |02,p). 
We note that the operator A can be viewed as one out- 
come of a generalized quantum measurement, albeit one 
that reduces the dimension of the measured system from 
four spins to one. The probabilities for these states are 
given by the decomposition (|14p , conditional on perform- 
ing this measurement and successfully obtaining the out- 
come associated with A. Fig. [1] illustrates this sampling. 
For sufficiently small pe, percolation theory tells us that 
the lattice decomposes into disconnected clusters with en- 
tanglement only within a cluster, and not between clus- 
ters. Specifically, if Pe is smaller than the critical bond 
percolation probability pj?°j"'^ for the lattice, then we ex- 
pect all the clusters will be sufficiently small to store and 
simulate efficiently on a classical computer. 



B. Classical resource requirements 



We now provide an explicit bound on the classical com- 
puting resources required to execute our simulation algo- 
rithm. Specifically, we prove that if Pe < pj?";"^, then the 
total simulation cost is bounded by a polynomial in N. 

In our simulation, each run yields a lattice with some 
entangled and some product-state bonds. We identify 
regions - "clusters" - that are connected by entangled 
bonds. Distinct clusters are therefore in a product state, 
and can be stored and simulated separately. For cluster 
Cj of size \Cj\, a direct simulation would have a cost pro- 
portional to 2^1"-^^ I . Because there are at most N clusters, 
the total cost of classically simulating a measurement se- 
quence on the state p{f3,9) is bounded above by N2'^^'~'^^ 
for each round of the simulation. 

We expect that for sufficiently high temperatures pe 
will become very small and the resulting clusters will 
be small enough to simulate efficiently. To confirm 
this expectation, we invoke results from percolation the- 
ory [1^, [l3| to bound \Cj\. The essential idea is that, 
if Pe is smaller than the critical bond percolation prob- 
ability Pc°it'^ for the lattice, then all the clusters will be 
"small enough" to simulate efficiently. Specifically, if 
Pe < Pcrit"^' then the mean cluster size x{Pe) docs not 
depend on A^. In addition, the largest cluster size is al- 
most surely of size ©(logjA^) with standard deviation 
©(logj log2 N) [13I. Now imagine we reserve N classical 
registers, each of kN'^ classical bits for some constants k 
and c, and use each register to store the state of one of the 
M (< N) distinct clusters. This allows us to store, with k 
bit precision, the quantum state of any cluster for which 
\Cj\ < log2 N'^, and to simulate local measurements on it 
efficiently. As the largest cluster is almost surely of this 
size, the total simulation cost is bounded by a polyno- 
mial in A^. Note that while this analysis shows efficient 
classical simulatability is possible, it is certainly not op- 
timal - more precise estimates would require considering 
the distribution of cluster sizes. 
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FIG. 2: The parameter space of our model, with bounds on 
the usefulness of the thermal state on a simple cubic lattice for 
MBQC. Region C represents the region on which any MBQC 
scheme performed on the thermal state can be efficiently simu- 
lated on a classical computer. Region Q represents the region 
where we can rigorously show, via the dephasing channel ar- 
gument given in the text, that the state is universal for MBQC 
using ideal measurements. Region Q' is a region where the 
local filtering yields a percolated thermal cluster state with 
T < 0.28A; it is possible to perform universal MBQC on any 
ground state or 9=0 thermal state in this phase, and states 
in this phase with T,9 > may be universal for MBQC al- 
though rigorous error correction thresholds for such states are 
not currently known. The black dashed curve represents the 
separability criteria of [31; above this curve, the PEPS de- 
scription of the thermal state is separable along any bipartite 
division of the cubic lattice given by a plane. 



improved lower bounds of this critical temperature can 
be obtained by allowing the value of uje to vary from bond 
to bond. 



IV. REGION Q: UNIVERSAL FOR MBQC 



C. Critical temperature for region C 

With this analysis, we can now lower bound the crit- 
ical temperature above which our simulation algorithm 
succeeds. In the case of 6* = 0, we can first compute the 
critical value of cUe = ^crit, and relate this to a critical 
temperature via 



tanh(/?eritA/2) = ]^t 



UJr 



(16) 



For some well known lattices we find that kTcru is 0.813A 
(Honeycomb), 1.6921A (Square), 7.1617A (Triangular), 
13. 1 A (Cubic). For 6' > 0, we resort to numerical meth- 
ods (see Appendix [B|) and again find a critical temper- 
ature. As expected, this temperature becomes zero at 
6 = 'n 12, when the ground state is a product state. Fig- 
ure [2] plots this solution for a cubic lattice. We note that 



We now demonstrate that, for the 'cold and weak field' 
region in parameter space, the thermal state of IH) is a 
universal resource for MBQC on an appropriate lattice. 
Specifically, we use the error thresholds of p, 0] j together 
with a local filtering method, to prove the existence of a 
finite region of parameter space for which the thermal 
state of ([T]) is a universal resource for MBQC given ideal 
single-qubit measurements. 

Consider performing the local measurement on every 
site in L as introduced in [lO| for the zero temperature 
case, described by the measurement operators 



Afo = \/l-tan2 0|O)(O|, hh = tan0|O)(O| + |1)(1| , 

(17) 
where 2(/) = 7r/2 — Q. For T = 0, the effect of the mea- 
surement is easily calculated by recalling that this state 
can be expressed as (ULc)\S)®''' ■ Both measurement op- 



erators Mq and Mi commute with CZc , and thus we can 
consider their effect on the state \6) at each site. If the 
'0' outcome is obtained, which occurs with probabihty 
Po = sin0, the post-measurement state is |0); if '1' is ob- 
tained, the resulting state is |-f-). If the operation CZc 
is apphed to a lattice of qubits where a subset C (cor- 
responding to those qubits for which the measurement 
outcome '1' was obtained) has each qubit prepared in 
the |-|~) state and all remaining qubits prepared in |0), 
the result is an ideal cluster state on C where the re- 
maining sites remain unentangled. We will refer to this 
procedure as local filtering. 

Similar results hold for the thermal state, with 
cos 2^ = tanh(/3A/2) smO. The resulting state is a ther- 
mal state of the 9 = Hamiltonian on C (the subset of 
the lattice where the measurement results '1' were ob- 
tained) with an increased temperature T' ~ {kBP')~^, 
given by 



tanh(/3'A/2) 



cos 9 



^l-tanh2(/3A/2)sin2e 



tanh(/3A/2). (18) 



The probability that, for a given site, the measurement 
yields '1' is pi = 1 - tanh(/3A/2) sin6l. 

With this filtering, we can argue the existence of a fi- 
nite region of parameter space where MBQC is definitely 
possible on a cubic lattice using ideal singlc-qubit mea- 
surements. Define C to be the subset of qubits for which 
the '0' outcome occurred, i.e., the complement of £'. For 
6* <C 1 the measurement disentangles the qubits in C 
from the rest of the cluster state, leaving them in the 
state |0). We then randomly flip each qubit in £', such 
that they are described by the completely mixed state. 
We also apply a Z gate to all neighbouring qubits to those 
that have been flipped into the state |1). The effect of 
this further processing is to prepare a = thermal clus- 
ter state for which complete dephasing has been applied 
qubits in C. We then discard the measurement record 
and the result is a thermal cluster state on the entire 
lattice C for which each qubit has been passed through 
an effective dephasing channel x{p) == (1 ^ p)p + pZpZ, 
with p = tanh(/3A/2) sin0. This is just the same state 
as the 9 = thermal state with T chosen such that 
P=(1+cxp(A/A:bT))-i. 

Thus, a thermal 6* 7^ state can be converted by local 
measurements into an ideal cluster state subjected to de- 
phasing noise. With this fact, we use the results of p], 
which demonstrate that MBQC can proceed using the 
cluster state on a body-centred cubic (BCC) lattice with 
dephasing noise up to pc ~ 2.9 x 10~^. This bound deflncs 
a region at low T and 9 wherein every state is useful for 
MBQC given ideal measurements. Fig. [5] shows this re- 
gion for a cubic lattice (noting that a cubic lattice can be 
converted into a BCC lattice with local Z measurements, 
which are unaffected by the noise), labeled as Region Q. 
At 9 = 0, the boundary corresponds to a temperature of 
Tc ~ 0.28A. 

However, these results may be too conservative. Con- 
sider the 9 y^ ground state, and apply the local fll- 



tcring measurement. Again, the effect of this measure- 
ment is to disentangle some qubits from the lattice, and 
leave the remaining qubits in a cluster state. One can 
then investigate whether the resulting subset C contains 
a cluster of neighbouring sites that spans the lattice £. 
This will occur with certainty if the success probability 
1 — sin6' is above the site percolation threshold, p^''°. So 
if l~sm9 > p^''°i then the resulting cluster state on C 
is a universal resource for MBQC [la, 13 • This lower 
bound (valid only at T = 0) is much larger than the con- 
servative lower bound on the critical value of 9 obtained 
via the dephasing argument above. 

A similar question can be asked of the thermal clus- 
ter state; however, the results of [3, M 13 do not directly 
apply because it is not clear how to convert the thermal 
cluster state on the irregular set C into the BCC lattice 
for which error thresholds are known. (A direct conver- 
sion using the methods of |16l . [l9[ would require perform- 
ing local X measurements, which do not commute with 
the noise.) As a result it is not clear that MBQC can 
in fact be performed. Nevertheless, Fig. [2] shows the re- 
gion Q' for which the success probability of the filtering 
operation is greater than the site percolation threshold 
1— sin 9 > Pc''°j with a resulting temperature T' less than 
the critical temperature for universal MBQC. An alter- 
native approach would be to apply the error correction 
procedure of [20| which corrects the effects of both finite 
temperature and qubit loss errors. The performance of 
this scheme in a three dimensional cluster is currently un- 
known, however it may yield better performance than the 
strategies described above in some regions of parameter 
space. 



V. CONCLUSIONS 

We have shown that the thermal equilibrium states 
of Eq. ID) undergo a transition in their usefulness for 
MBQC, from a region of parameter space where every 
state is a universal resource, to one where every state 
is efficiently classically simulatable. In spite of this dra- 
matic change in computational power, these states do not 
exhibit any corresponding phase transition. 

We conclude by contrasting our results with some re- 
lated work. First, in a similar model with a transverse 
rather than longitudinal local field, the system does ex- 
hibit a zero-temperature phase transition at sufficiently 
high field strength, and correlation functions can be iden- 
tified which characterize a phase for which a set of MBQC 
gates can be performed over arbitrary ranges [2ll |. Thus, 
a transition in computational power coincides with a 
quantum phase transition in this model. 

Second, non-analytic behaviour in long-range entan- 
glement quantities such as localizable entanglement [l| 
does not necessarily indicate a phase transition. In our 
model, for sin 6' < 1— p^'*'', the localizable entanglement 
in the ground state is precisely equal to 1 ebit at all length 
scales, and at 6* = 7r/2 it is equal to zero (because the 



state is a product state). There is no analytic function 
fulfilling these requirements, and so the localizable entan- 
glement is a non-analytic function of 6, indicating a sharp 
transition somewhere in the interval 1— p^'*'' < sin^ < 1. 
Again, we emphasize that there is no quantum phase 
transition at any value of 9 in this model, and so we have 
an example of a system where a non-analyticity in the 
localizable entanglement length does not identify an un- 
derlying quantum phase transition. This in contrast to 
the situation for a large number of models discussed in 
the literature [l| in cluding the cluster Hamiltonian with 
a local X field [2l|,[2l. 
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the projection. Consequently, we cannot sample directly 
from the ensemble of pure bond states according to the 
probabilities in the decomposition of Eq. p^ : instead, we 
use a sampling procedure where each step is conditioned 
on the success of the site projections. Here we describe an 
algorithm that samples efficiently from the appropriate 
distribution, and show below that the algorithm samples 
the correct distribution for the whole system. The key 
requirement for this simulation to be efficient is that the 
success of the measurement Ai at site i is independent of 
the other sites in the lattice. 

The following algorithm efficiently samples from the 
posterior bond distribution. That is, it reproduces the 
distribution p ({p^*-'^ } \Ai (g) A2 (8) . . . (8) An) of a par- 
ticular configuration {p^^^^} of all bonds in the lat- 
tice conditional on success of each of the projectors Ai. 



set all virtual qubits to empty 
for all bonds {i,j) do 

set any empty qubits in {i. * \i.j} and {j. * \j.i} to 

the state po = tr2 [pbond] 

sample from the distribution 

p[p^'''^\A, ® A,,p^"\'-^\p^'-*\^-'^] 



APPENDIX A: THE SAMPLING ALGORITHM 

The sampling of PEPS bonds to generate a single 
pure state instance of the thermal state requires some 
care. We must sample from the ensemble of bond in- 
stances after each of the site projections A = |0) (0 . . . 0| + 
|1)(1 . . . 1|. To obtain the correct distribution, the pro- 
jections A should be viewed as physical measurement op- 
erations, and one should sample from the posterior dis- 
tribution conditional on the success of the operations. 
However, the success of the projection at a site can vary 
depending on the choice of pure bond states at the site. 
For example, a site containing a virtual qubit in the state 
|0), and another virtual qubit in |1), will yield zero con- 
ditional probability for the projection, while a site with 
all qubits in the pure state |0) will guarantee success for 

I 



set the state of qubits i.j and j.i to the corresponding 
bond state 
end for 

Here, {i. * \i-j} is the set of virtual qubits at site i 
except for the one associated with the bond {i,j), and 
p[pi'-i'i\Ai(»Aj,p^'-*\'-^\pi^-*\^-'^] is the posterior prob- 
ability for the bond {i,j), given the state of the other 
virtual qubits at sites i and j and given that the projec- 
tions Ai and Aj succeed. We can calculate this distribu- 
tion straightforwardly from Baycs' rule. 

We now show that the above algorithm indeed samples 
from the correct probability distribution, which is the 
posterior distribution for the bond configuration, condi- 
tioned on the success of all N projections. The desired 
distribution may be written, using Bayes' rule, as 



P 



({/'^">}l^i 



•A2 



, .«) Ajv 



p{Ai®A2...®An\ {p^*^> }) X p ({p{*^> }) 



p (Ai (g) ^2 «) • ■ • «> An) 
On the other hand, the algorithm described above samples, on a bond-by-bond basis, from the distribution 



(AI) 



Palg 



{i.*\i.j} p{j-*\j-i}] 



[pi^^^}) = l[p[p^^^^A(^A,,i 
{y} 

p[A, ® Aj\p^'^\ p^"\'-^\ p^^-*\^-'^ X p{*J'> 



n 



p[Ai (g) Aj\pi^-*\'-i},pii-*\^-^}] 



(A2) 



Here, the product runs over all bonds {ij} in the lattice, p[Ai (g) Aj\p^^^\ pU-*V-J} ^ p{j*\j*}] jg the probability that the 
projections at site i and j succeed, conditioned on the bond {ij} being in state p^*-'^ , and p^*-'^ is the prior distribution 
for each bond, as given by the probabilities {pe,Ps,fj,} of Eq. p4|) . The denominator p[Ai (g) ^j|p^*-*\'-^^,p^-'*\^*^] 
denotes the total probability that the projections at site i and j succeed, conditioned only on the state of the virtual 
qubits {i. * \i.j} and {j. * \j-i}- 

In order to show the equivalence of the two distributions of Eq. (JA1[) and Eq. (|A2p , first consider the denominators 
in Eq. (|A2p . These may be written as 

p[A, (g A,\p^''-*\'-'\p^^-*\'-''^] = ^p[A, ® A,>i'^'\p{"V-j}^pb-.*b-}] X pi''^ , (A3) 

k 

where pj,^ arc the elements of the ensemble decomposition of Eq. ()14p . and p^*"' the corresponding weights. This 
expression may be written 

p[A, (g) Aj|p{"V-^'>, p{J-*\^-^>] = ^ tr[(A, <g ^,>{"V-''> (g pI''^ <g p{J-*\^-^>(A, (g ^j)^]pi'^'^ 

fc 

= tr[(A, ® A,)p{"\^-^'> ® p^^2d ® P^^-*\^'-'^(^. ® Aj)^] , (A4) 

where p^ond denotes the total state of the bond. In Sec. [El we show that for 9 > 0, the thermal bond state is of the 
form pbond = \{I + ct^ ^ Z + jZ (g /)(/ + aZ g) X + 7/ g) Z), where a and 7 are parameters, determined numerically 
for a given temperature and field strength (see Sec. |B]). For any physical values of a, 7, p{'*\*-'l and p{J-*\-J'l, one 
finds that this expression may be written 

p[A, ® A,|p{"V-^},pb-*\j-.0] ^ p[A,|p{"\-^">,pJ'-^>] X p[A,\pi^-*\^-'\pl'-'^] . (A5) 

The numerators in Eq. (|A2p contain likelihood factors of the form 

p[A, (g Aj\p^'^\ p^'-*\'-^\ p^^-*\^-'^ = tr[(A, g) AjOp^''*^*'^'^ «) P^*^'^ ® p^'-*\^-'^A^ g) y4j)^] . (A6) 

If p^*-'^ = p^'-'^ eg p^^-^^ then it is clear that this expression can be written as a product of two probabili- 
ties (corresponding to independent outcomes Ai and Aj at sites i and j respectively). Conversely, as we show 
in the final section of this appendix, the entangled component of the decomposition is also of the form pg = 
jil + ctoX g) Z + -foZ g) /)(/ + uqZ (g at + 70/ g) Z), where ao and 70 are real parameters. Since, in this case, 

pbj} is of the same form as the state p^ond appearing in Eq. (jA4p . it follows from Eq. ljASP that the expression can 
also be written as a product of two probabilities. Thus this term may be written as the product: 

p[A^ ® Aj\p^'^\pi"\'-^\p^J-*\^-'^] = p[A,|p{'-J>,p-f'-*\*-^'>] X p[Aj\p^J-'\pi^-*\^-'^] . (A7) 

Note that generalizations of Eq. (|A5[) and Eq. (|A7p to arbitrary numbers of sites also hold. For example, in the 
three site case, where site j neighbours both site i and site fc, we find 

p[A, (g Aj (g Afe|p-f'-*\*-^'>,p{^-*\{J-'J-'^'»,p{'=-*\'=-^'>] 

= p[A|p{"V-^>,p^;'-^>] X p[A,\p^=-*\^^-^'^-'^^\pi=-'\pl,'-'^] X p[A,|p{'=-A'=-J},p{'=-^}] , (A8) 

and also 

p[A,g)Aj®Afe|p{'^>,p{j''^>,p{"*V-j}^p{j-*\b-.j-fc}}^^{fc.*\fc.j}] 

= p[A,|p{*-^'>,p{'-*V'-J}] X p[Aj|p{J-*>,p{^-'^>,p{j-*\{^-''^-'^»] X p[Afe|p{'^-^>,p{'^--*\'^--J>] . (A9) 



One way to verify these expressions, and their generalizations to any number of sites, is as follows. Recall that the 
site projections of the PEPS representation commute with CZc gates applied on each of the bonds. Consider a single 
instance of the thermal ensemble of entangled and product state bonds, as shown in Figure [T]Jb). Consider performing 
a CZc operation on every bond that is occupied by an entangled state, leaving the bonds with product states alone. 
The site projections are unaffected by this operation, however the entangled bond states have become completely 
unentangled and the resulting state is a product state between all sites. The success probabilities at each site are now 
clearly independent. 



Palg 



Using the factorizcd expressions Eq. (jASp and Eq. (jA7[) . Eq. (|A2p becomes 
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p[A,|p{-^^>,pp>, . . . ,p^'-^^>] X p[A,|p{^-^^},p{'-^^},pp> 



.P^^^lx 



xp[A,|p{'-Ji>,p{--'"^},...,p{-J'<^>] 



p[A,|p^'-^^>,p^^-^"^>, . . . ,p^;'-^^>] X p[A,\pi^-n},p'i;-^^\. . . ,p{'-^''>] X ... X p[A,\p{^-n}^p{^.n}^ . . . ,p{'-^^}] 
{*^> (AlO) 



xJIp 



where in the second Hne we have re-arranged the product such that each factor now corresponds to a particular site of 
the lattice, rather than a bond. The ordering of the virtual qubits at site i implied by the indices ji . . . jd corresponds 
to the order in which the corresponding bonds are sampled in the algorithm. Eq. (jAlOp can be simplified considerably 
by noting that many repeated terms appear in both the numerator and denominator, yielding 
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Finally, by observing that the expres- 
sions p \Ai ® A2...® An\ Ip^'-'H) and 
p {Ai ® A2® . . .® Am) of Eq. (|A1|) may be written in 
a factorized form (by making use of the appropriate 
generalization of Eqs. (jASp and (jA9[) . we can identify 
each term in Eq. (jAlip with the corresponding term in 
Eq. (JA1[) . Thus, our algorithm indeed samples from the 
correct distribution for the whole state. 



APPENDIX B: THERMAL STATES WITH 
NON-ZERO LOCAL FIELD 

In the non-zero field case the thermal bonds are deter- 
mined by the requirement that they project, under Ai at 
each site, to the correct thermal cluster state. 

In the zero temperature case it is possible to obtain an 
analytical expression for this bond state. It is straight- 
forward to see that the states 

Phond{T = 0,6') 

= - (/ + a: ® Z cos 20 - Z ® / sin 2(/)) 

X (/-hZ(g)Xcos20-/(g)Zsin2(/)), (Bl) 

yield the correct pure cluster state, provided that 
tan'^(0 4- 7r/4) — tan(0/2 -I- 7r/4), where d is the coordi- 
nation number of the lattice. However, for finite temper- 
atures it is slightly more difficult to obtain the relations 
between the bond parameters and the values of T and 



9. As stated in the main text, the bond states take the 
general form 

Pbond(T,6') = -{I+aX®Z+-iZ®I){I+aZ®X+-^I®Z) , 

(B2) 
where the parameters a and 7 obey 



a = tanh 



/3A 



cos^ 2_] 



J even 



j odd ^-^^ 



/5A . 



tanh -— sin 9 \^ 



J even 



(B3) 



while also being constrained to producing physical bond 
states. In one dimension these equations are relatively 
straightforward and may be solved without much trou- 
ble, however for the cubic lattice the conditions involve 
sixth order polynomials in 7 and so we opt to solve 
these constraints numerically, and find that for any T 
or < ^ < 7r/2 we obtain an appropriate thermal bond 
state. 

To calculate the classical simulation bound from per- 
colation on the lattice, the thermal bond state must be 
decomposed into an entangled part and a separable part. 
Instead of calculating the Best Separable Approxima- 
tion (BSA) for the thermal state, we use the state p^ = 
Pbond(2^ = 0,0) as the (non-maximally) entangled pure 
bond state in the decomposition pbond = PePe + (1 —Pe)Ps 
and choose the largest Pe for which p^ remains separable. 
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